Correlated cluster mean-field theory for spin systems 
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A cluster mean-field method is introduced and the applications to the Ising and Heisenberg models 
are demonstrated. We divide the lattice sites into clusters whose size and shape are selected so that 
the equivalence of all sites in a cluster is preserved. Since the strength of interactions of a cluster with 
its surrounding clusters is strongly dependent on the spin configuration of the central cluster itself, 
we include this contribution in the effective fields acting on the spins. The effects of "correlations" 
between clusters can be taken into account beyond the standard mean-field level, and as a result 
our cluster-based method gives qualitatively (and even quantitatively) correct results for the both 
Ising and Ifeisenberg models. Especially, for the Ising model on the honeycomb and square lattices, 
the calculated results of the critical temperature are very close (overestimated by only less than 5 
%) to the exact values. 



PACS numbers: 75.10.-b, 05.50.+q 



I. INTRODUCTION 

In general, it is difficult to find the exact solution of 
an interacting many-body system, except for a few sim- 
ple cases. Hence, various mean-field theoriesi^i^i^i^ have 
been proposed and widely used to get some insight into 
the properties and behavior of the systems. A mean- 
field concept was first used by Weissi for the problem of 
ferromagnetism. Weiss's molecular-field theory (MFT) 
was successful in providing a qualitative understanding of 
the transition from ferromagnetic to paramagnetic states. 
After the initial success by Weissj^i the MFT has played 
an important role in the studies of many-body systems 
and is still often used as a starting point of the investi- 
gation due to its simplicity. 

In the vicinity of the critical point, mean- field- type 
theories generally breaks down since statistical fluctua- 
tions neglected in them become important. Therefore, in 
order to study the detailed critical properties (e.g., crit- 
ical exponents), it is necessary to resort to other tech- 
niques, such as the renormalization-group methodsi^'^'^ 
However, since mean- field- type methods are capable of 
predicting the behavior of the system over a wide range 
of parameters (e.g., temperature) at a relatively low cost, 
a further development of those methods is still impor- 
tant, and many attempts have been made even in recent 
years^iii^ii^ 

The Ising model is the simplest nontrivial model of 
cooperative phenomena. The Hamiltonian for this model 
is given by 



(1) 



where the sum runs over all nearest-neighbor pairs, J is 
the coupling strength, and the spin Ui takes either -fl 
or —1. In the MFT, the interactions between a spin and 
its surrounding spins are treated approximately, and the 
many-body problem is reduced to an effective one-body 
problem of a single spin in the effective magnetic field 
/leff = zJm, where z is the number of nearest neighbors 



and m is the magnetization of the system. Solving the 
self-consistent equation for the magnetization m — (a), 
where (• • • ) is an average over the ensemble, one can ob- 
tain the well-known result for the critical temperature of 
the transition from ferromagnetic to paramagnetic states: 



ksTc/J = z. 



(2) 



Mean-field-type methods generally give only classical 
predictions for the critical exponents, but the value of 
the calculated critical temperature is a reasonable mea- 
sure for assessing the accuracy of the method. 

To obtain more accurate results than this "one-site 
cluster" approximation, it is quite natural that one at- 
tempts to increase the size of clusters i^'^i^ii^'i^ii^'i^iii 
For example, in the so-called Bethe-Peierls- Weiss (BPW) 



approximation 
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using a cluster of (z -\- 1) sites (one 



central spin and its surrounding z spins), one can par- 
tially take into account the effects of spin correlations and 
fluctuations. The interactions of the central spin with its 
nearest neighbors (called the "first shell" ) are treated ex- 
actly, while the influence of spins outside the cluster is 
replaced by an effective field hes which acts on the z spins 
of the first shell. Unlike the MFT, the effective field hcg 
is determined by the condition that the average value of 
the central spin should be equal to that of the spin on a 
first-shell site. Recently, Du et al}^ extended the BPW 
method by using a group of chains composed of a central 
chain and its nearest-neighbor chains instead of the clus- 
ter of {z -\- 1) sites. Also, Etxebarria et alJ^ proposed 
another extended BPW method and demonstrated that 
it yields a fairly accurate estimate of critical temperature 
for the square-lattice Ising model. 

Oguchi's method^ii^ is a more straightforward way to 
improve the results of the MFT. In this method, one con- 
siders a cluster consisting of Nc neighboring spins. The 
interactions between the spins at the cluster edge and 
the outside spins are treated approximately as effective 
internal fields. Calculating the average (ai) for each site 
in the cluster, one regards the mean value of them as 
the magnetization of the system, i.e., m — (l/Nc) '^(o'i). 
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Then the effective fields are determined self-consistently 
by the value of m. Clearly, the case of iVj, = 1 corresponds 
to the conventional MFT. Although it is expected that 
one can obtain closer results to the exact values as the 
size of the cluster Nc is increased, this method has two 
shortcomings. First, one has to deal with a quite large 
cluster to get sufficiently accurate results. In practice, it 
is difficult to solve exactly the A^c-site problem of such 
a large cluster. Second, since the system is divided into 
finite-size clusters, equivalence of all sites (or periodicity 
of the system) is generally lost with few exceptions. 

Recently, Wysin and Kaplar*^ made a significant im- 
provement to the MFT in a quite simple way. Their 
"self-consistent" correlated field (SCCF) approximation 
remains basically the same as the standard MFT in the 
sense that one deals with effective one-site problems. 
In this method, the effective field acting on a central 
spin from its surrounding spins takes two different values 
h^g = zJm^, depending on the state of the spin itself 
(cT; = +1 or —1). Here m+ {m^} is the average value of 
the spin on a first-shell site when the state of the cen- 
tral spin is fixed to be -1-1 (—1). Taking into account 
the effects of correlations between neighboring spins in 
this way, they obtained more accurate results of criti- 
cal points compared with some other methods, such as 
the BPW approximation and the Onsager reaction field 
(ORF) correction.aaoas 

In this paper, as an attempt to further improve these 
methods, we introduce a new cluster-based approxima- 
tion method, which we refer to as the "correlated cluster 
mean-field" (CCMF) theory. Depending on the lattice 
type, we select a cluster of a different size and shape. 
For example, the calculation for a square lattice model is 
based on square-shaped four-site (2 x 2) clusters. Clearly, 
the four sites in this cluster are equivalent by symme- 
try. Additionally, to take into account the effects of the 
cluster-cluster correlations beyond the standard mean- 
field level, we use a similar idea as in the SCCF method 
of Ref . S 

The paper is organized as follows. First, in Sec. [TTl 
we demonstrate the application of our modified cluster 
mean- field theory to the Ising model on several lattices. 
In Sec. lIIIi the accuracy of our method is verified by com- 
paring the obtained results with those obtained by some 
other methods. In Sec. IIVI the cluster-size dependence 
of the results is discussed. Then, in Sec. |Vl taking the 
isotropic Heisenberg model in a uniform field as an exam- 
ple, we extend our approach to quantum spin systems. 
Finally, a summary is presented in Sec. IVII 



II. CORRELATED CLUSTER MEAN-FIELD 
THEORY FOR THE ISING MODEL 

In the following, we shall demonstrate the application 
of our method to the Ising model on typical four types 
of lattices (honeycomb, triangular, square, and simple 
cubic) one by one. 




FIG. 1: Schematic representation of the honeycomb lattice 
divided into hexagonal clusters. The arrows indicate the ef- 
fective fields acting on the sites (a) in cluster C and (b) in 
cluster C' when the spin at site 1 is fixed to be -1-1 or —1. 



A. Honeycomb lattice 

When the coordination number z is 3, unfortunately, 
the SCCF approximation yields a wrong result: the equa- 
tion for Tc has no solutioni^ It is therefore important to 
verify whether the critical point can be calculated by ap- 
plying our method to the case of the two-dimensional 
(2D) honeycomb lattice. We first divide the lattice into 
six-site hexagonal clusters as seen in Fig. [Ija). The 
strength of interactions of a cluster with its surrounding 
clusters should be strongly dependent on the spin config- 
uration of the central cluster itself. Thus, including this 
contribution into the effective fields acting on the spins, 
we consider the following six-site problem on cluster C: 

Hc^-J ^^^J ~ (3) 
(ij)GC iec 

where the first (second) sum runs over all nearest- 
neighbor pairs (all sites) within cluster C, and h"^^ = 
Jm'^' is the effective field acting on site i from the neigh- 
boring spin in the nearby connected cluster (namely, e.g., 
from the spin at site 4' for i = 1). Here m'^' is the mean 
field of the neighboring spin of site i and we assume that 
its value depends on the state of spin z, 

= ( S'^' = (4) 
[ m (ctj = -1). ^ ' 

Especially in this case, since there is only one bond be- 
tween two clusters, the situation is similar to that in the 
SCCF approximation. Note that the values of m+ and 
m~ do not depend on the site number i due to the sym- 
metry of the hexagonal cluster. 

When the values of to+ and m~ are obtained, the six- 
site problem described in Eq. ([3]) is exactly solvable, and 
we can calculate the magnetization given by 

m= (d,) =Tr(a,e-^^^)/Tr(e-'5^^), (ieC), (5) 
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FIG. 2; Schematic representation of the square lattice divided 
into square-shaped clusters. The arrows indicate the effective 
fields acting on the sites (a) in cluster C and (b) in cluster C' 
when the spins at site 1 and site 2 are fixed to be +1 or —1, 
respectively. 
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(8) 



It is clear from the symmetry that the effective field from 
site 3' to site 2 should be h^^"^ = Jm"''-''^ [see Fig. dja) 
for the others]. In this way, within our approximation, 
the states of the four spins in cluster C can be described 
by the Hamiltonian 



He = -J o-icr 



"-off '^i 



3 

J (o-i(T2 + 0-20-3 + 0'30'4 + O'40'i) 



+ rn° 



icri 

I (72 

1 0-3 



(9) 



where /3 = l/ksT. 

The values of m"*" and m~ are determined in the follow- 
ing self-consistent way. Let us now focus on the cluster 
next to C, which we call cluster C" [see Fig. (Ub)]. One 
obtains the values of m+ and by calculating the aver- 
age values of the spin at, for example, site 4' in Fig.^b) 
when the spin at site 1 is fixed to be and —1, respec- 
tively, i.e.. 



(^4')|,^^±i =Tr(a4 



/Tr 



where 



-J 



E 



(6) 



(7) 



The upper (lower) signs correspond to the case where the 
value of cTi is fixed to be -1-1 (—1). Thus it is only neces- 
sary to solve the set of two equations (namely, two six-site 
problems) given by Eqs. ([6]) and ([7]) self-consistently for 
TO"*" and m~ and then one can calculate the magnetiza- 
tion m from Eq. 



B. Square lattice 

In the case of the square lattice, there are two bonds 
between two neighboring clusters as shown in Fig. [2] and 
this makes matters somewhat more complicated. The 
spins in cluster C" are strongly affected by the states 
of the spins at site 1 and site 2 of cluster C, and thus, 
including the effect of the feedback, the effective fields 
acting on site 1 and site 2 should depend on the states 
of themselves. Then, for example, the effective field from 
site 4' to site 1 can be denoted by h'^^'^'^ = Jmf^'^'^ , where 



where X^igc denotes the sum over all nearest neighbors 
of site i within cluster C . The values of the four mean 
fields shown in Eq. ([5]) can be determined by 



\ ^ ' \ a\ —S,(72—S 



— Tr ( e 



(a4'e-^^c/)/Tr(e-''^c/) , (10) 



where 



jCTi (77 

fteff Oi 



({M}#{3',4'},{4',3'}) 

= —J (cri'0'2' -f CF2'(7y + 0'3'(T4' -I- cr4'(Ti') 
-J(m'^i""2' +m"^"'^')av 

-J{m'^^"^3' ) (72' 

-J(m'"3"^2' +s')cr3, 

-J(m'"4""i' +s)cr4'. (11) 



The notation H^f, corresponds to the case where the val- 
ues of (Ti and a2 are fixed to be s and s', respectively 
(s,s' — +1 or —1). Solving the set of four equations 
given by Eqs. pl}|) and pT|) self-consistently, one can ob- 



tain the values of the mean fields 



", and 



m . Then the magnetization m [Eq. ([5])] can be calcu- 
lated by solving the four-site problem given by Eq. ([5]). 



C. Triangular lattice 

For the triangular lattice, we perform the calcula- 
tion based on triangle-shaped three-site clusters shown 
in Fig. [3] in a basically similar fashion to the previous 
two cases. From the geometry, there are two types of 
effective fields acting on a site. Let us specifically focus 
on site 1 in cluster C. There are four bonds from sites 
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FIG. 3: Schematic representation of the triangular lattice di- 
vided into triangle-shaped clusters, (a) The effective fields 
acting on site 1. (b) Positional relation between three clus- 
ters (C, C", and C"). 



in other clusters to site 1: two of them are from clus- 
ter C" and the other two are from cluster C" and C'-'^-', 
respectively. 

Since cluster C is connected to cluster C at only one 
point, the spins at site 2' and site 3' are directly affected 
only by the state of the spin at site 1. Thus we assume 
that the two effective fields from sites 2' and 3' depends 
only on the state of the spin at site 1: 2/iJ^j^ = 2Jm°'i 
[see, Eq. ^ for the definition of m'^']. 

As for cluster C" [C^'^^), on the other hand, there are 
connections with cluster C at two points, site 1 and site 
3 (site 1 and site 2). Thus, in our procedure the effective 
field from cluster C" {C'^^^) to site 1 should depend on 
the spins at these two sites: /i^" = Jm"^''^ (/iX' = 
Jm'^^"^) [see, Eq. ^ for the definition of m'^''^^]. When 
taken together, the total effective field acting on site 1 is 
given by 2hl^ + h^p + h^' [see Fig. EJa)]. 

Eventually, for the case of the triangular lattice, we ob- 
tain the following Hamiltonian describing the three spins 
in cluster C: 

He = ~J ^'"^J - H^/ioff^'' - J2 KP""^ 

= - J ((Tio-2 + cr20-3 + o-30-i) 

-J (2m'"i + m"^^"" + nf^''^)ai 
-J {2rrf^ + nf'"^ + nf^''^)a2 
-J {im"^ + m'"-""i + m'"-""^) erg. (12) 

The magnetization m can be calculated from this Hamil- 
tonian via Eq. 

In a similar way to the previous two cases, let us deter- 
mine the values of the mean fields m+, mT ^ m++, , 
m '", and m . First, fixing the value of ui and consid- 
ering the three-site problem of cluster C", one can obtain 
the equations for the mean fields and m~ , 

-^ = (-2')L,^±i=Tr(a2'e-'^«c') /Tr(e-^^c^) ,(13) 



where 

- Y Ka'^i T J (o-2' + C^S') 

({»J}/{2',3'},{3',2'}) 

= — J (cri'cr2' + f2'0'3' + Ca'CTi') 

-J (2m'"2' -f to'^2"^i' ± 1) 0-2/ 

-J(2m'"3' +to''-^""i' ± 1)0-3.. (14) 

The upper (lower) signs correspond to the case where the 
value of (Ti is fixed to be -t-1 (—1). Due to the symme- 
try, of course, one can obtain the same value of by 
calculating (0-3') |^^^_|_^ instead of Eq. 

Next, we calculate the other type of the mean fields: 
771"*"+, m"' , m '", and m . To this end we now focus 
on, for example, cluster C" [see Fig. [H^b)]. In a similar 
manner to the case of the square lattice, we can derive 
the following equations for the four mean fields: 



= Tr (a^^.e-f'^c'') /Tr (g-'^^c'') , (15) 

where 

Hc" ^ -J Y ^^'^3 ^ X! ^Ks^^ 

{i.3)ec" tec" 

(»#2") 

- ^ Ks^'^i - J {3 + s') 0-2" 

i,ieC" 

= — J ((Ti"0'2" + f2"(T3" + Cr3"(Ti") 

-J(2to'"i" +771"^""^-" +m'"i"'^3")cri// 

-J {m"^""^" + m"^""^" + s + s') (72" 

-J {27X1"^" -f m'"3"'^i" +m'"3"'^2")cr3„. (16) 

In the case of the triangular lattice, one have to solve 
the set of six nonlinear equations in terms of the six pa- 
rameters m+, TO~, m^"*", vi^ , m '", and m , which 
are given by Eqs. P^ - fTB]) . Nonetheless, since the clus- 
ter we selected here consists of only three sites, it is easy 
to solve the set of equations. 

D. Simple cubic lattice 

As an example of three-dimensional (3D) lattices, we 
show the application of our method to the case of the sim- 
ple cubic lattice, where the calculation is performed based 
on eight-site cubic clusters. In this case, clusters are con- 
nected by four points of each face (see Fig. |4]), and thus 
we have 2^ mean fields [m'^^'^'^ ^rri^'^'^ , . . . ,m ). 
For example, the effective field from site 1' to site 2 is 
defined by a function of the states of the four spins cr2. 
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TABLE I: Reduced critical temperatures {kBTc/J) for the Ising model on various types of lattices from different approximations 
and exact or series values. 



Lattice 


MFT(z) 


Exact or series 


BPW 


ORE 


SCCE 


SME 


CCME 


Honeycomb 


3 


1.519 


1.820 






1.464 


1.593 


Square 


4 


2.269 


2.885 




2.595 


2.142 


2.362 


Triangular 


6 


3.641 


4.933 




4.788 


2.628 (3.543") 


4.519 


Simple cubic 


6 


4.510 


4.933 


3.955 


4.788 


4.570 


4.753 



"From the corrected SMF with a = 0.8 (see Ref.Hal). 



TABLE IL Comparison of the results obtained by the CCMF method for the square-lattice Ising model with those obtained 
by the BPW method, its extended version proposed by Etxebarria et al. (Ref. 11^ ). Kikuchi's square approximation (Refs. 
and [2^) ■ and exact values. Ci{Tc) and C2(Tc) are the values of the correlation functions between first and second neighbors, 
respectively, at Tc- Nc is the size of the cluster, and M is the number of nonlinear equations involved. 

Method M kgTc/J CijTc) C-zin) 

Exact 2.269 0.707 0.637 

BPW 5 1 2.885 0.333 0.111 

Etxebarria et al. 12 4 2.351 0.607 0.501 

Kikuchi 4 3 2.426 0.562 0.438 

CCMF 4 4 2.362 0.608 0.495 
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FIG. 4: Schematic representation of eight-site clusters for the 
simple cubic lattice. 



(T3, (Tg, and (T7. Due to symmetries, only 12 of the 2^ 
mean fields are actually independent. The calculation is 
straightforward extension of that in the case of the square 
lattice, so we display only the final result in Sec. IIIII 



III. DISCUSSION OF RESULTS FOR THE 
ISING MODEL 

Before presenting the results obtained by the above 
procedures, we will briefly review the screened magnetic 
fleld (SMF) approximation, which was recently proposed 
by Zhuravlevii^ In his method, one considers an effective 
one-site problem of a certain spin (the central spin) in an 
effective magnetic field like the standard MFT. The influ- 
ence of the interactions between the central spin and all 
others should decrease quickly with the distance between 
them. In the SMF, he took into account this "screen- 
ing effect" in the following way. Introducing the factor 

2/2 

/o-o^ he assumed that the effective field acting on the 



central site can be represented by 

/leff = Jm^e-'-'/'^o, (17) 

where the sum is over all spins except for the central spin 
itself. As a characteristic length oq, the lattice constant 
was chosen in Ref. [ij. For the calculation of the critical 
points, one have only to replace the effective field in the 
MFT by Eq. Then one obtains 

M/J-^e-'-'/"". (18) 

In Table HI we summarize the results of the critical 
temperature Tc obtained by the CCMF and other meth- 
ods. We can see that the results obtained by our CCMF 
method are very close to the exact results or approximate 
value from the series-expansion method;^'^i2^ Partic- 
ularly, for the honeycomb and square lattices, our re- 
sults {kBTc/J = 1.593 and 2.362) overestimate the cor- 
responding exact values by only less than 5 %. For the 
triangular lattice, the accuracy is inferior as compared 
with the above two cases, but is better than those of the 
other methods [except for the "corrected" SMF (Ref. ol) 
with the correction coefficient a — 0.8]. 

In most previous approaches, such as MFT, BPW, and 
SCCF, the results are dependent only on the coordination 
number z. In contrast, our cluster approach can distin- 
guish between two lattices with different geometries but 
equal coordination numbers, as well as the SMF method. 
For example, the obtained critical points of the triangu- 
lar and simple-cubic lattices are different, although they 
have the same coordination number z = 6. This is be- 
cause we selected different clusters depending on the ge- 
ometries of the lattices. For the 3D simple cubic lattice, 
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unfortunately, our result of the critical point represents 
only a little improvement from that of the SCCF ap- 
proximation. This is attributed to the fact that the im- 
portance of including cluster correlations is relatively low 
compared with 2D cases, since mean-field approaches are 
generally speaking expected to work better in higher di- 
mensions. 

Next, focusing on the case of the square lattice, let 
us compare the results obtained by the CCMF method 
with those obtained by the extended version of the BPW 
method proposed by Etxebarria et alr^ and Kikuchi's 
square approximation.i^*^ In the method of Etxebarria 
et al., a relatively large cluster (consisting of 12 spins) 
is considered, and one introduces one effective field and 
three effective couplings between the "boundary" spins, 
which are determined by the condition that the period- 
icity of the system is preserved (for detail, see Ref. Il2h . 
The so-called cluster-variation method (CVM) proposed 
by Kikuchi* provides a systematic hierarchy of approx- 
imations, each of which is specified by the set of "basic 
clusters" used for the calculation. For example, in the 
"square" approximation^'-'^ of the CVM, the largest ba- 
sic cluster is a 2 X 2 square. 

In Table [Tll we summarize the results of the critical 
temperature Tc and values of the connected correlation 
functions between first and second neighbors at Tc, ob- 
tained by four different cluster-based approaches. In 
the CCMF calculation, the correlation functions between 
site i and site j can be obtained by {(JiCfj) — {cFi)"^ ~ 



Tr {a,Uj exp[-/37?c]) /Tr (exp[-/3i/c]) 



(*,J e C). 



The extended BPW approximation gives the most accu- 
rate value {UbTc/ J = 2.351) compared with the other 
approximation methods in Tables [J and |TT1 although one 
have to deal with a 12-site problem. Kikuchi's square ap- 
proximation gives UbTc/ J = 2.426. The CVM results are 
usually more and more accurate as the size of the consid- 
ered basic clusters increases. Using the Tanoji approxi- 
mation,^'' where a nine-site cluster is used as the largest 
basic cluster and 44 variational parameters (M — 44) 
are needed, one can obtain the more accurate result 
{kBTc/J = 2.346). Our CCMF procedure can achieve 
the same level of accuracy by treating only four sites. 
(However, note that in the CCMF approach it is neces- 
sary to solve a set of four slightly different four-site prob- 
lems to obtain the values of all unknown mean fields.) 



IV. CLUSTER-SIZE DEPENDENCE AND 
CRITICAL EXPONENTS 

In Sec. ini we have demonstrated the application of 
the CCMF method to the Ising model by using a cluster 
of appropriate size and shape for each lattice type. For 
the square lattice, we employed the clusters illustrated in 
Fig. O Hereafter we refer to this procedure as the square 
approximation of the CCMF method for descriptive pur- 
poses. 

In Table IIIII we show the values of fc^Tc/J for the 



TABLE III: The values of ksTc/J and the mean-field criti- 
cal coefficients for tfie square-lattice Ising model in the single 
site (SCCF), pair, and square approximations of the CCMF 
method. 



Approximation ksTd J 


nis Jx+ Jx~ 


SCCF 2.595225 
Pair 2.437857 
Square 2.362099 


2.534519 0.716621 0.358311 
2.922637 1.130168 0.565084 
3.551663 1.537636 0.768818 


TABLE IV: Comparison of the CAM estimates and exact val- 
ues of critical exponents. 




P 7 = 7' 


CAM 
(fcsT.VJ) 


0.343 2.138 
(2.331) (2.118) 


Exact 2.269 


0.125 1.75 



square-lattice Ising model in three different approxima- 
tions: the single site (i.e., SCCF), pair, and square ap- 
proximations of the CCMF method. In the "pair" ap- 
proximation, we adopt a nearest-neighbor pair of spins 
as a basic cluster for the calculation (for detail, see the 
Appendix). It can be seen that the value of the crit- 
ical temperature becomes more and more close to the 
exact one (fc^Tc/J = 2.269185 ••■) with increasing the 
size of the used cluster. However, in all the cases, the 
obtained critical exponents are classical ones, as in most 
other mean-field approaches. Namely, the temperature 
dependence of the spontaneous magnetization just below 
Tc is given by 



T -Tr 



Tr. 



1/2 



for T <Tr. 



(19) 



and the zero-field susceptibility just above T^ and just 
below Tc are given by 



X- 



T-Tr 
Tr 
T-Tr 



Tr 



for T>Tr 



for T <Tr. 



(20) 



The values of the critical coefficients m and x± are also 
shown in Table IIIII 

The coherent-anomaly method (CAM) proposed by 
Suzuki^^ allows one to estimate non-classical critical ex- 
ponents from a series of mean-field-type approximations. 
Now, assuming that the approximations in the CCMF 
scheme constitute a part of a canonical series^ let us es- 
timate the non-classical critical exponents for the square- 
lattice Ising model. According to the CAM, the critical 
coefficients to and x± should behave asymptotically as 



m ^ ci<5(r,)-(V2-/3), 
X- ^ c,SiTrr^'''-'\ 



(21a) 
(21b) 
(21c) 



7 



for 



0, 



(22) 



where {ci} are certain constants, T* denotes the exact 
value of the critical temperature, and Tc is the approx- 
imate critical temperature obtained in each approxima- 
tion. Since each equation has three unknown variables, 
Ci, T* , and the critical exponent {f3, 7, or 7'), we need the 
approximate values of Tc and the corresponding critical 
coefficients (to, x+, or x~) for, at least, three different 
levels of approximations. 

In Table |TVl we show the values of P, 7, and 7' ob- 
tained from Eq. (PT|) and the data obtained from the 
three levels of approximations of the CCMF method pre- 
sented in Table Hm Although non-classical exponents are 
obtained, the values of them are not very accurate. This 
may be because the size of clusters used in the approx- 
imations is relatively small, It is expected that more 
accurate results are obtained by taking a larger cluster 
like that used in the Tanoji approximation of the CVM^ 
However, since the initial lattice symmetry is artificially 
broken, i^'^^^ the CCMF theory is not directly applicable 
for a cluster larger than a certain critical size (e.g., the 
2x2 cluster for the square lattice) in the present form. 



APPLICATION TO THE HEISENBERG 
MODEL 



In this section, we consider the application of the 
CCMF theory to a spin-i Heisenberg ferromagnet in a 
uniform magnetic field as an example of quantum spin 
systems. The Hamiltonian of the system is given by 



(23) 



Here Si = {Sf , Sf , Sf) is the usual spin operator at site 
i, which satisfies the commutation relation [S'f,S'j'] = 
iifiuxS^Sij, = Sf ± iSf are the spin raising and low- 
ering operators, and h denotes an applied magnetic field. 

One of the advantage of cluster-based approaches, such 
as the BPW and CCMF methods, is that the effects of the 
spin-flip term, Sf Sj +S^Sj', can be taken into account 
by the straightforward application. Now, we focus again 
on the case of the square lattice, and calculate the magne- 
tization of the system by using the CCMF approach. In 
the square approximation (see Fig. [J), we first introduce 




0.0 



0.5 1.0 1.5 

temperature ksT/J 



FIG. 5: The temperature dependence of the magnetization 
of an isotropic Heisenberg ferromagnet at h/ J — 0.1. Com- 
parison of the result obtained from the CCMF approach with 
those of the MFT, BPW, RPA, Kikuchi's square approxima- 
tion, and the QMC calculations. The error bar of the QMC 
results are smaller than the line width. 



the four effective fields h. 



s 



where 



1 o- 

1 
2 
1 
2 
1 



s: 



s: 



1 

2+^ 
1 
2 
1 
2 
1 



5*^ 



2 



m 



{(7i CTj =T), 

{(7i =T, (Tj =i), 
(CTi =i, (Tj =i). 



(24) 



Here, ui =] , [ represents the z component of the spin at 
site i. The procedure for the calculations of these mean 
fields and the magnetization to is almost the same as that 
in Sec. IIIBI The difference is only that the Ising inter- 
action between neighboring spins in a cluster is replaced 
by the Heisenberg one. Si ■ Sj, in this case. 

The temperature dependence of the magnetization 
TO obtained by the CCMF approach is compared with 
those from the MFT, BPW, the so-called random- 
phase approximation (RPA) j^^i^^ Kikuchi's square 
approximation,— and quantum Monte Carlo (QMC) cal- 
culations^i2£i^i2^ in Fig. [5l The BPW approxima- 
tion, which was first applied to the isotropic Heisenberg 
model by Weiss in Ref. [l^, gives a qualitatively incor- 
rect behavior in the low-temperature region. This be- 
havior becomes more pronounced at lower magnetic field 
strengths. Therefore, when applied to other systems, the 
BPW method has a risk that one cannot decide whether 
the behavior of the obtained results is an essential prop- 
erty of the system or an artifact of the approximation. 
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The similar behavior in the low temperature region is 
also observed in the result of Kikuchi's approximation. 
In the "Bethe-type" or "effective-field" approaches, such 
as the BPW method^ i^^'^^ and its extended versions, ^^'■^^ 
the values of the effective fields (and effective couplings) 
are determined by the periodic boundary conditions and 
thus these approaches generally overestimate the corre- 
lation effects between the spins. 

The Green's function method with the RPA decou- 
pling approximation (or "Tyablikov decoupling" method) 
is one of the most popular approximation schemes to deal 
with quantum spin systems. As is shown in the figure, 
the RPA is a good approximation over a wide temper- 
ature range. However, since the treatment for the Ising 
term of the Hamiltonian is identical to that of the MPT, 
the approximation should get worse as the Ising-type 
anisotropy of the system increases. Especially, when the 
RPA scheme is applied to the Ising model, the obtained 
results are the same as that of the MFT. Besides, a spe- 
cial attention is required in applying the RPA scheme to 
a relatively complicated system as is pointed out in our 
recent study. '^'^ 

In the "Weiss-type" approaches, the effective fields 
are defined in terms of some kind of mean fields (e.g., 
hcs = zJni in the simple MFT), and those values are 
determined by the corresponding self-consistency con- 
ditions. Oguchi's method, the SCCF, and CCMF ap- 
proaches are categorized into this group. These ap- 
proaches generally underestimate the correlation effects 
between the spins and therefore overestimate the magne- 
tization. The CCMF method also overestimates the mag- 
netization as is shown in Fig. [51 but the difference from 
the QMC data is much smaller than that of the MFT. 
Moreover, above a certain temperature (fcsT/J ~ 0.7), 
the accuracy exceeds that of the RPA scheme. The 
results obtained here and in Sec. IIIII indicate that our 
CCMF method can give qualitatively (and even quanti- 
tatively) correct results over a wide range of temperature 
and strength of the exchange anisotropy. 



VI. SUMMARY 

In this paper, we have proposed a cluster-based mean- 
field approach, which we refer to as the CCMF theory 
and have demonstrated its applications to the Ising and 
isotropic Heisenberg models as examples. By using the 
CCMF approach, one can take into account the effects 
of spin correlations between clusters beyond the stan- 
dard mean-field level. Since a cluster of different size 
and shape is used for each lattice type, the obtained re- 
sults are dependent not only on the coordination number 
z but also on the geometry of the lattice. 

The results obtained from the CCMF method are in 
good agreement with the corresponding exact values, se- 
ries expansion, or QMC data for the both Ising and 
isotropic Heisenberg models. Especially, for the Ising 
model on the honeycomb and square lattices, the cal- 



culated results of the critical temperature are very close 
(overestimated by only a few percent) to the exact ones. 

As well as being very accurate, the advantage of this 
method is that it is widely applicable. Since the cal- 
culation in the CCMF approach is based on clusters 
consisting of several spins, in addition to the Ising-type 
interaction, the contributions of other types of interac- 
tions (e.g., the Heisenberg-type exchange interaction as 
we demonstrated) can be taken into account in a straight- 
forward way. Although we have considered here only 
the two simple cases, this method is expected to be use- 
ful also for studies of more complicated and interesting 
systems. For example, the CCMF method should work 
well for the XXZ model (Heisenberg-Ising model) with 
easy-plane anisotropy^SiM and, of course, with Ising-type 
anisotropy, which lies between the two models we consid- 
ered here. Also, the extensions to systems with four-spin 
(ring) exchange interactions,'^- higher spins, and random 
systems^SiSl should be considered in future studies. 
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APPENDIX A: PAIR APPROXIMATION OF 
CCMF 

In the pair approximation for the square-lattice Ising 
model, the states of a nearest-neighbor pair of spins (see 
Fig. [H]) are described by the Hamiltonian 

He = -J J2 ''i'^j-^Kk^i- Yl '^KP<^^ 

-J {m"^ +2771"^"^) (72- (Al) 

Oguchi's pair approximation"" is obtained by requiring 
the equality to"" = m'^'^ = m, where m — (1/2) X^c^'''')' 
in the above Hamiltonian. In the CCMF method, the 
values of the mean fields are determined by solving the 
set of self-consistent equations 

m± = (a2')L^=±i=Tr(a2'e-'='^c') /Tt {e~^"c'^ (A2) 

and 



= Tr (ar.e^'^^c'') /Tr (^e-^"c''^ , (A3) 
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FIG. 6: Schematic representation of the square lattice divided 
into pairs of neighboring sites. The arrows indicate the effec- 
tive fields acting on the sites in cluster C. 



H. 



(ij)ec" leC" i,iQC 

—sJayf — s' Ja2" 
= -Jay, (72" - J [m"^" + ni'^^'"'^" + s) ay 
-Jim"^" +m'"2"'^i" +s')(72". (A5) 



where 



= — J(Ti'(T2' ~ J {ra 1 + 2m i ^ ) ui' Then the magnetization m is calculated from Eqs. 

-J(2m'"2"'i' ± l)cr2' (A4) and jJIl)- 
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